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Abstract 

The BCH formula of Rezek and Kosloff is a convenient tool to 
handle a family of density matrices, which occurs in the study of 
quantum heat engines. We prove the formula using a known argument 
from Lie theory. 

1 Introduction 

The Hamiltonian 

H = -hu(a^a + aa^) (1) 

of the quantum harmonic oscillator belongs to the Lie algebra su(l,l) with 
generators 

51 = |((« f ) 2 + a 2 ), 

5 2 = ^((a f ) 2 -a 2 ), 

5 3 = -(a^a + aa^). (2) 

As a consequence, it is possible to write down simplified Baker-Campbell- 
Haussdorf (BCH) relations [UI2]- These have been used to study the quantum 
harmonic oscillator in a time-dependent external field [21 IU El El CO El E]- 
The topic of the present paper is a new BCH relation, introduced recently 
by Rezek and Kosloff [10]. They consider the family of density matrices 

P "^(/3,7) () 



1 



with real /3 and complex 7, and with 

Z{p,i)= Tre 7a2 e-^e 7(at)2 . 



(4) 



Because all operators appearing in fl3]) belong to the Lie algebra it is clear 
from the general Baker-Campbell-Haussdorf relation that it must be possible 
to write 

The explicit expression of the coefficients \ and £ as a function of j3 and 7 
is found in the Appendix of [10]. The functions were derived [11] using the 
algebraic manipulation software Mathematica. 

Note that the special case of fl5J) with £ = appeared in the physics 
literature before (see Example I of Section II of [T2]; see also [131 El])- In the 
present paper the relation fl5]) is derived using the argument of [12] . 

The relation (pj]) is of interest in its own. But it is also very useful in the 
study of quasi-stationary processes JTUJ [15] . Indeed, from (JSJ) it is clear that 
the density matrix p describes a system in thermal equilibrium at inverse 
temperature £. On the other hand, the expression (j3J) is more convenient for 
practical calculations. 

We derive the BCH formula in the next section. In Section 3 follows 
a similar BCH formula valid for su(2). In the final section follows a short 
discussion. 



2 The identity 

The the l.h.s. of (jSJ) can be written in terms of the generators of the Lie 
algebra su(l,l) as 

e 2 7 (5i+iS , 2) e -2 ( 9^53 e 27(Si-iS'2)_ (g) 

These generators satisfy the commutation relations 

[Sx,S 2 ] = 1S3, 
[£2, £3] = -iSi, 

[S 3 ,S 1 ] = -iS 2 . (7) 

Introduce su(2) generators 7\ = —iSi, T 2 = iS 2 , and T 3 = S3. Then (jHJ) 
becomes 

X = e 2 7(m+T 2 ) e -2/3r t a;r3 e 27(iT 1 -T2) > /g\ 
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The relation (jSj) does not depend on the choice of the representation of the 
su(2) algebra. Therefore, we may change it. A favourable choice is that of 
the Pauli spin matrices o a = 2T a . Using that {a\ ± ia 2 ) 2 = and a 2 = I the 
calculation becomes very easy. One obtains 



X — gj7(ci— 10-2) ^-fiftwoz gi7(cri+i<T2) 

= (I + ij(ai — icr 2 ))(cosh(/j,a;) — a 3 sinh(/zu;))(I + vf((Ti + 2cr 2 )) 
= e"^ 3 - 2k| 7 | 2 + zk( 7 + t>! + k( 7 - t> 2 + 2k| 7 | V 3 , 



with k, = e ^ as before. On the other hand is 



(9) 



exp [xa 2 - + x(a f ) 2 ) 



exp (2x(Si + ^>S 2 ) - 2f fajS 3 + 2*(Si - ifi^)) 
exp (2x(iTi + T 2 ) - 2£HujT 3 + 2x{iT l - T 2 )) . 

(10) 



In the Pauli spin representation this becomes e Y with 
Y = i(x + X)°i + (x ~ X)o2 ~ Zfwa 3 . 



Because the Pauli matrices anti-commute and their squares equal I there 
follows that 



Y 2 = A 2 I with A = ^(M 2 — 4|x| 2 - 
Hence one obtains 

e Y = cosh(A) + - sinh(A)F. 
A 

Comparison with OH]) gives the 4 conditions 



cosh(A) 
^sinh(A)(x + x) 



A 



sinh(A)(x - x) 



a + k, 
«(7 + 7) 9 
«(7-7), 



A 



sinh(A)£ftu; = a, 



(12) 

(13) 

(14) 
(15) 

(16) 

(17) 



with 



a = smh(/3hu) — 2k\j\ 



2k 



[l - k 2 -4k| 7 | 2 ] . 



;is) 



The solution of these equations is 



no; smh(A) 

X = «-4r7. (20) 
smh(A) 



with 



sinh(A) = ^a 2 -4k 2 | 7 | 2 (21) 



These results coincide with those found in the Appendix of [10J. 

Note that the expressions for £ and x can be inverted easily. Given £ and 
X one obtains A from (j!2p . Then a follows by inverting (1191) . This gives 

. sinh(A) . . 

a = Hu£ -^A (22) 

A 

Next P is obtained from (|14|) 

k = cosh(A) — a. (23) 

Finally, 7 follows from (T2"U|) 

sinh(A) 

7 = Px- 24 

3 An example with SU(2) symmetry 

Formulas similar to (J5J) can be derived for other symmetry groups than 
SU(1,1). For instance, in the case of SU(2) one has 



e 



i° +e -^et°- = exp ( xa+ - & z + x<r_) (25) 



with a± = \(<J X ± icr y ). Using a\ = 0, a\ = I, <j±<t z = T°± ; an d 0+0"- = 
|(1 + cr 2 ) the l.h.s. becomes 

l.h.s. = (1 + 7cr+) (cosh(/3) - sinh(/?)a 2 ) (1 + 7<r_) 

= cosh(/3) + -|t| V + ^(70+ + jaJ) - (sinh(/3) - -e p \l\ 2 )a z . 

(26) 

The r.h.s. of fl25|) is evaluated using 

( X a + -^ + x^-) 2 = A 2 I, (27) 
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with A = a/£ 2 + Ixl 2 - One finds 

r.h.s. = cosh(A) + — sinh(A) (x a + ~ i a z + X a ) ■ 

A 

Equating both expressions yields the set of equations 



(28) 




cosh(A), 



— — sinh(A)£ 
A 

- sinh(A)x- 



(29) 



(30) 




(31) 



Given £ and %, the value of A can be obtained from its definition. The 
solution then reads 



Conversely, given f3 and 7 one obtains A from fT29|) . Then £ and x follow from 
([30]) and fl3Tl) . respectively 

4 Discussion 

The BCH relation of Rezek and Kosloff is somewhat special because it is 
written in a form suited for application to density matrices. Similar results 

found in the literature [TJ[2j[3j|lJ[5j[6j[71|8j|9] a * m a ^ ^ ne ca l cu l a ^ on °f 
time evolution operators and refer to similarity transformations, this is, to 
expressions of the form e A Be~ A . But the l.h.s. of (J5]) is not a similarity 
transformation. This is precisely the reason why this BCH relation is of 
interest! The change of the spectrum implies that the average energy (H) = 
Tr pH will depend on the value of the parameter 7. This dependence is 
essential in the context of heat engines. 
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